On the nature of steady states of spin distributions in the presence of spin-orbit 

interactions 
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In the presence of spin-orbit interactions, the steady state established for spin distributions in an 
electric field is qualitatively different from the steady state for charge distributions. This is primarily 
because the steady state established for spin distributions involves spin precession due to spin-orbit 
coupling. We demonstrate in this work that the spin density matrix in an external electric field 
acquires two corrections with different dependencies on the characteristic momentum scattering 
time. One part is associated with conserved spins, diverges in the clean limit and is responsible 
for the establishment of a steady-state spin density in electric fields. Another part is associated 
with precessing spins, is finite in the clean limit and is responsible for the establishment of spin 
currents in electric fields. Scattering between these distributions has important consequences for 
spin dynamics and spin-related effects in general, and explains some recent puzzling observations, 
which are captured by our unified theory. 
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I. INTRODUCTION 

Spin-orbit interactions are frequently the most impor- 
tant factor determining spin dynamics in solids. From a 
technological perspective, novel physical phenomena that 
may lead to improved memory devices and advances in 
quantum information processing have been shown to be 
intimately related to spin-orbit interactions^ Practical 
applications usually rely on generating and maintaining 
a spin polarization, and in this context spin-orbit 
interactions can play a constructive or a destructive role. 
In an external electric field spin precession gives rise to 
steady-state spin densities^^&Zi 8 .' 9 - ' 10 ! 11 ! 12 and spin 
currents ; 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 ' 23 - 24 - 25 - 26 - 27 - 28 - 29 - 30 - 31 - 32 - 33 



quency fife /h = [fife [/ft and is scattered to a different 
wave vector within a characteristic momentum scatter- 
ing time T p . Throughout this paper we assume that 
£FT p /h ^> 1, where Ef is the Fermi energy, which is tan- 
tamount to assuming that the carriers' mean free path is 
much larger that their de Broglie wavelength. Within 
this range, the relative magnitude of the spin preces- 
sion frequency fi& and inverse scattering time l/r p de- 
fine three qualitatively different regimes. In the ballis- 
tic (clean) regime no scattering occurs and the temper- 
ature tends to absolute zero, so that eft p — ^ oo and 
QkTp/h — ► oo. The weak scattering regime is character- 
ized by fast spin precession and little momentum scatter- 
34 , iftgedffis to.9 ami S'slisM.mimms&ks .temperature , yielding 



yet spin precession is often the leading cause of spin po- 
larization decay i 51 i 52 i 53 i 54 These facts suggest that spin 
precession plays a nontrivial role in the establishment of 
a steady state for spin distributions in electric fields. 

Spin-orbit interactions are present in the band struc- 
ture and in potentials due to impurity distributions. 
Spin-orbit coupling is in principle always present in im- 
purity potentials and gives rise to skew scattering, which 
leads to the extrinsic anomalous Hall effect and the ex- 
trinsic spin-Hall effect^ Band structure spin-orbit cou- 
pling may arise from the inversion asymmetry of the 
underlying crystal lattice^ (bulk inversion asymmetry), 
from the inversion asymmetry of the confining poten- 
tial in two dimensions^ (structure inversion asymme- 
try), and may be present also in inversion symmetric 
systems.^ 8 Band structure spin-orbit coupling gives rise 
to spin precession, is the cause of magnetic anisotropy 
in magnetic semiconductors^ and causes spin flips even 
in the course of elastic scattering by spin-independent 
potentialsj 53 i 60 Band structure spin-orbit interactions in 
spin- 1/2 electron systems can always be represented by 
a Zeeman-like Hamiltonian H = (1/2) er • flk describ- 
ing the interaction of the spin with an effective wave 
vector-dependent magnetic field flk- An electron spin 
at wave vector k precesses about this field with fre- 



EFTp/h » flkTp/h 3> 1. In the strong momentum scat- 
tering regime £FT p /h > 1 > Oji t p /H. 

In this paper, we will be concerned with the inter- 
play of external fields, scattering processes, and band 
structure spin-orbit interactions in establishing a steady 
state for spin distributions. We will concentrate on 
electrons in uniform electric fields. In charge transport, 
the steady state is characterized by a nonequilibrium 
correction to the density matrix that is divergent in 
the clean limit, indicating a competition between the 
electric field, accelerating charge carriers, and scattering, 
which inhibits their forward motion. On the other hand, 
nonequilibrium corrections that arise as a result of band 
structure spin-orbit coupling in crystal Hamiltonians re- 
present a different kind of interplay between the electric 
field and scattering processes. Firstly, the spin-orbit 
splitting of the bands gives rise to spin-dependent scat- 
tering even from spin-independent potentials. Secondly, 
the presence of spin precession causes the steady state 
established for spin distributions to be highly nontrivial. 
We will demonstrate that the density matrix contains 
a contribution due to precessing spins and one due to 
conserved spins, and the steady states established for 
each of these are qualitatively different. Steady state 
corrections oc t p , which diverge in the clean limit, are 
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associated with the absence of spin precession and give 
rise to spin densities in external helds 1 1 1 1 1 1 2 

Steady state corrections independent of t p , which 
are finite in the clean limit, are associated with spin 
precession and give rise to spin currents in external 
fields . 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28 - 29 - 30 - 31 - 32 - 33 - 34 



tributions of conserved and precessing spins, and attempt 
to shed light on the physical content of vertex corrections 
in systems with spin-orbit interactions. In the process, 
we give a detailed exposition of the underlying theory, 
which was also employed in our recent worki^ 
S^^^^M 2 4M^-^§^M^k follows. In the next 



Scattering between these two distributions induces sig- 
nificant corrections to steady-state spin densities and 
spin currents. 

In research on electrical generation of spin den- 
sities and currents two rather different approaches 
are employed. Linear response theories based on 
Green's functions^^^ 4 .^^!^^ provide a dia- 
grammatic interpretation of spin-related phenomena in 
electric fields. Such theories usually provide reliable re- 
sults, but the physical picture is not always evident. 
Theories based on a kinetic equation for the density 
matri^^ i 39 ' 40 ! 41 ' 42 tend to be transparent. However, 
it is often difficult, in such theories, to determine from 
the outset all the terms that play an important role in 
the processes under study. Correspondences between the 
two approaches mentioned were identified in a recent pa- 
per by Sinitsyn et al^- The formalism used in the present 
paper, although relying on a density matrix, is equivalent 
to Kubo linear response theory. 

We consider large, uniform systems, working in mo- 
mentum space without making semiclassical approxima- 
tions (by which we understand approximations pertain- 
ing to simultaneous consideration of a particle's position 
and momentum.) The role of scattering in spin-related 
effects is an issue that is currently under intense inves- 
tigation. In our general formalism, the terms responsi- 
ble for scattering are derived rigorously beginning with 
the quantum Liouville equation for the density operator. 
One particular issue at the center of current debate is 
the physical interpretation of vertex corrections in spin- 
related phenomena. In charge transport, it is well known 
that vertex corrections decrease the weight of small-angle 
scattering. We demonstrate that, in spin transport, scat- 
tering phenomena have two main consequences. Firstly, 
the driving term in the equation for the spin density ac- 
quires a contribution due to spin-dependent scattering. 
Secondly, scattering between the distributions of con- 
served and precessing spins enters the equations for these 
two distributions. Our work suggests that both of these 
processes are contained in vertex corrections. 

In a recent related article^ we focused on the spin 
current response of a semiconductor to an electric field, 
showing that spin currents are not restricted to the spin- 
Hall effect and dwelling briefly on the distinction between 
nonequilibrium spin currents and spin densities excited 
by electric fields. In this article, we concentrate on the 
steady state and on the difference between the steady 
state established for spin distributions and that estab- 
lished for charge distributions. We investigate aspects 
of the steady state for spin distributions that arise as a 
result of spin precession and have no analog in charge 
distributions. We provide a detailed analysis of the dis- 



section we start from the quantum Liouville equation 
and derive an equation describing the time evolution of 
the electron density matrix including the full scattering 
term in the first Born approximation. Section III is de- 
voted to the complex case of steady states in the presence 
of band structure spin-orbit coupling, demonstrating the 
existence of two spin distributions and their subtle inter- 
play. Analytical expressions are given for general scat- 
tering as far as possible, as well as exact solutions for 
short-range impurity scattering. We close with a sum- 
mary of our findings. 



II. TIME EVOLUTION OF THE DENSITY 
OPERATOR 

We consider a system of non-interacting spin-1/2 elec- 
trons. The electrons are represented by a one-particle 
density operator p. The expectation value of an observ- 
able represented by a Hermitian operator O is given by 
tr(pO), which motivates us to study the density opera- 
tor p in detail. The dynamics of p are described by the 
quantum Liouville equation, which is projected onto a 
set of states of definite wave vector (in which the matrix 
elements of p form the density matrix) , and an equation 
is obtained for the time evolution of the density matrix. 
This equation is valid for any elastic scattering in the first 
Born approximation. In this article we discuss systems 
with long mean free paths and do not consider diffusion 
terms explicitly. Consequently, semiclassical approxima- 
tions are not necessary, and the method employed in the 
present work is equivalent to the Kubo linear response 
formalism. 



A. Quantum Liouville equation 



The quantum Liouville equation satisfied by p is 



dp i ~ - 

^ + - [H + u,p] = a. 



(i) 



The Hamiltonian H, considered here in the framework 
of the envelope function approximation (EFA), contains 
contributions due to the kinetic energy and spin-orbit 
coupling. The effect of the lattice-periodic potential of 
the ions is taken into account through a replacement of 
the carrier mass by the effective mass. We shall hence- 
forth refer to H as the EFA Hamiltonian. Scattering 
is introduced into the system through the potential U, 
which may be due to impurities, phonons, surface rough- 
ness, or other perturbations. In this article we focus on 
impurity scattering, as the effects that are discussed are 
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frequently observed at very low temperatures, where the 
role of phonon scattering may be neglected. 

The Liouville equation is projected onto a set of time- 
independent states of definite wave vector {\ks)}, which 
are not assumed to be eigenstates of the EFA Hamilto- 
nian H. The matrix elements of p in this basis will be 
written as pkk' = Pkk' = (ks\p\k's'), with correspond- 
ing notations for the matrix elements of H and U. Spin 
indices will not be shown explicitly in our subsequent 
derivation, the quantities pkk' > Hkk' , and Ukk' being 
treated as matrices in spin space, pkk' is referred to as 
the density matrix. With our choice of basis functions of 
definite wave vector, matrix elements of the Hamiltonian 
Hkk' — Hk 5kk> are diagonal in k. However, if the EFA 
Hamiltonian contains spin-orbit coupling terms, the ma- 
trix elements Hk are generally off-diagonal in spin space. 
Matrix elements of the scattering potential Ukk' are off- 
diagonal in k . Matrix elements diagonal in k in the scat- 
tering potential would lead to a redefinition of Hk , which 
is analogous, in Green's function formalisms, to the offset 
introduced by the real part of the self energy. Scattering 
is assumed elastic and, given that we will work in the 
first Born approximation, the matrix Ukk' is assumed 
diagonal in spin. (We thus do not take into account 
so-called skew scattering terms, which require terms of 
third order in Ukk' as well as explicit inclusion of spin- 
orbit coupling in the scattering potential, and which were 
studied in a recent work42) As mentioned above, we fo- 
cus in this article on impurity scattering. The impurities 
are assumed uncorrelated and our normalization is such 
that (ks\U\k's')(k's'\U\ks) = ni\Ukk'\ 2 ^ss' , where rij is 
the impurity density. Ukk' refers therefore to the matrix 
elements of the potential due to a single impurity. Ex- 
plicit expressions for these matrix elements for a screened 
Coulomb potential in two and three dimensions are given 
in Appendix A. 

Pkk' is divided into a part diagonal in k and a part 
off-diagonal in fc, given by p k k' = fk <W + 9kk', where, 
in gkk', it is understood that k k'. The quantum 
Liouville equation can be broken down into equations for 
fk and g k k' 

^ + ^[H k Jk] = -^[U,g]kk, (2a) 
^T + ^ H ^]kk' = -- h [Uj + g] kk '- (2b) 



In the first Born approximation the solution to Eq. (|2bD 
can be written as 



gkk' 



dt'e 



-iHt' 



U, f(t - 1') 



AHt' | 

' \kk'- 



(3) 



In order to shorten the equations, factors of h that ap- 
pear in the time evolution operators will be omitted, i.e., 



AHt' 



AHt' j% 



These factors will be restored in the fi- 



nal results. Since EFT p /h 3> 1, we shall expand f(t — t') 
in the time integral around t and, noting that terms be- 
yond f(t) are of higher order in the scattering potential, 
we shall only retain the first term, fit). The equation 
for fk then becomes 



^ + \[H k Jk] + J{fk)=Q, 



in which the scattering term J(fk) is given by 



(4a) 



1 f°° 



U,e 



-iHt' 



u,f(t) 



iHt' 



kk 



(4b) 



The integral over time in Eq. (|4b[) can be performed by 
inserting a regularizing factor e~ nt and letting 77 — > 
subsequently. We remark that, for potentials diagonal in 
spin space and spin-degenerate bands, Eq. (|4b[) simpli- 
fies to the customary expression for Fermi's golden rule. 
Therefore, Eq. Q can be viewed as a generalization of 
Fermi's golden rule that explicitly takes into account the 
spin degree of freedom. The commutator present in Eq. 
(|4a[) is a commutator in spin space that represents the 
effect of spin precession due to spin-orbit interactions. 



B. Scattering term 

The scattering term J(fk) will now be evaluated for 
an electron system with spin-orbit interactions. After 
inserting a complete set of states, J{fk) becomes^ 



/>oo 

J(fk) = 3 hm / dt'e-^' V \U kk/ e-^' t '(Uk'kfk - fk'U k 'k) e tHkt ' - e~ lHkt ' (U kk ' fk' - fkUkk')e lHk>t ' U h , k • 

( 5 ) 

The EFA Hamiltonian contains a kinetic energy term and a spin-orbit coupling term, Hk — H^ lu + In spin- 1/2 
electron systems, band structure spin-orbit coupling can always be represented as a Zccman-likc interaction of the spin 
with a wave vector-dependent effective magnetic field thus = (1/2) er • Qk- Common examples of effective 
fields are the Rashba spin-orbit interactional which is often dominant in quantum wells with inversion asymmetry, 
and the Dresselhaus spin-orbit interaction, 5 - which is due to the inversion asymmetry of the underlying crystal lattice. 



4 



The spin-orbit interaction for electrons is usually much smaller than the kinetic energy at typical Fermi energies, with 
the result that terms that are second order in the ratio of the two, Hj^/H^ ln , can usually be ignored. In the basis in 
spin space spanned by spin eigenstates | f) and | j) (commonly referred to as the Pauli basis), taking into account 
the fact that Ukk' is diagonal in spin space, the scattering term simplifies to 



J(fk) 



h 2 



lim 



d d k' 
(2-K) d 



w kk , 



fk') 



-iH k t' 



(fk - fk')e 



iH k ,t' 



(6) 



Here Wkk' = |C^fcfc'| 2 is the transition rate, and sums over wave vector have been converted into integrals following 
the standard procedure J2k' ~* Vd f d d k' /(2n d ), where d is the dimensionality of the system and the normalization 
volume Vd is chosen to be the d-dimcnsional unit cell volume. The density matrix fk is decomposed into a scalar part 
and a spin-dependent part, fk = nk 1 + Sk ■ Performing the time integral, after a series of lengthy but straightforward 
calculations, which are summarized in Appendix B, the scattering term can be expressed in the form (Jo + J s ) (rik) + 
Jo(Sk), with 



Jo (Xk) 



■KtiiVd f d d k' 



2h 



J T0f Wkk ' {Xk ~ Xk,) ^ (efc+ ~ e ' k 



S(e k - 



_) + S(e k - - e' k+ )] , (7a) 



Js(Xk) = 



^ 7n — w Wkk' (Xk — Xk') c • (fifc 

(Z7T )" 



2h 
2h 



n k ,)[S(e k+ -e' k+ )-6(e k --e k _)] 



J jSf Wkk ' {Xk X k ' ) fx (O fc - ft,,,) [6(e k+ -4-)- 6(e k - - 4+)] 



(7b) 



Xk represents either nk or Sk, and flk is a unit vector 
along flk- The energies tk± represent the eigenvalues 
of Hk, and are given by €k± = £ofc ± where the 

kinetic energy eofc = H 2 k 2 /2m* and fl k = \il k \. The term 
J s (X k ) in Eq. ([7b|l illustrates the fact that, when spin- 
orbit interactions are present in the band structure, even 
a spin-independent scattering potential usually gives rise 
to spin-dependent scattering in the scattering integral. 



Time evolution of the density matrix in an 
external electric field 



In the presence of a constant uniform electric field E, 
k = q — eEt/h is the gauge-invariant crystal momentum 
(with q the canonical momentum.) The states \ks) are 
chosen to have the form \ks) = e lq ' r \uk s ), where \u ks ) are 
lattice-periodic functions. We subdivide fk = fok + fEk, 
where the equilibrium density matrix fok is given by the 
Fermi-Dirac distribution, and the correction fEk is due 
to the electric field E. To first order in E, the correction 
fEk satisfies 



Ek 



[H, f E k] + J (fEk) 



eE df 



Ok 



(8) 



dt h l ,J ^' w """ h dk 
The matrix fEk is in turn divided, as above, into a scalar 
part and a spin-dependent part, fEk = n-Bfcl + ^Bfe- 
To first order in H^/H k ln , the scattering term can be 

expressed as J (fEk) = (Jo + Js) (risk) + Jo (Ssk), where 
the scattering operators have been defined in Eq. {?])■ 

The driving term arising from the electric field in Eq. 
© is (eE/h) -dfok/dk. The equilibrium density matrix, 



/ofc, is subdivided as fo k = fiofcl + Soki with a corre- 
sponding subdivision for the driving term. The equation 
for n E k is 



dn E k , i , v eE dn ok 

-dr +Jo{nEk) = ir--dk- 



(9) 



The solution of this equation is given by the well-known 
expression 

eEr p dn ok , . 

UEk = h~ ' -flfc"' (10) 

in other words, nEk describes the shift of the Fermi 
sphere in the presence of the electric field E. The ex- 
pression for the momentum relaxation time t p is a little 
different depending on the dimensionality of the system. 
In three dimensions, using 7 to denote the relative angle 
between k and fc', 



1 mkv d=3 r 



r d=3 



2TTfl 3 



d^ sin7 Wkk' (1 — COS7). (Ha) 



In two dimensions, with the same notation for 7, 



1 



mV d = 



r d = 2 2ttH 3 J 



d-fWkk' (I-COS7). (lib) 



The spin-dependent part of the nonequilibrium correc- 
tion to the density matrix SEk is interpreted as the spin 
density induced by E. The equation governing the time 
evolution of SEk is 

dS E k . i ,„ c , * , c , eE dS ok * / v 

— r T Wk, JEk\ + Jo WEk) — "Z - 1 —qT Js [nEk)- 

(12) 
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It is seen from Eq. (|T2]) that spin-dependent scattering 
gives rise to a renormalization of the driving term in the 
equation for SEk- Evidently, this renormalization has no 
analog in charge transport, see Eq. ©. 

So far, our work has not been restricted to the steady 
state. In this context, we remark briefly that Eq. (fT2|) 
can be used to describe spin relaxation, and is valid 
both in the presence and in the absence of external elec- 
tric fields^ (In the presence of electric fields SEk is an 
electric-field-induced nonequilibrium correction, whereas 
in their absence it is to be interpreted more generally as 
a nonequilibrium correction.) 



III. STEADY STATES FOR CONSERVED AND 
PRECESSING SPINS 

In the presence of band structure spin-orbit interac- 
tions, an electron spin at wave vector k precesses about 
an effective magnetic field f2fc. The spin can be resolved 
into components parallel and perpendicular to Qk- In 
the course of spin precession the component of the spin 
parallel to is conserved, while the perpendicular com- 
ponent is continually changing. It will prove useful in our 
analysis to divide the spin distribution into a part repre- 
senting conserved spin and a part representing precessing 
spin. This is accomplished below. 



A. Distributions of conserved and precessing spins 

Firstly, the effective source term, which enters the RHS 
ofEq. (fT2|) . is divided into two parts, (eE/h)-dSok/dk~ 



J s (riEk) = ^Ek\\ Here, 'Sem commutes with the 

spin-orbit Hamiltonian and is given by T*Ek\\ — a fc^fc°j 



where 



tr 



eE 



Ok 



Ok 



J s in E k) 



Hi 



(13) 



while T.Ek± is the remainder. In matrix language Y,Ek± 
is orthogonal to the spin-orbit Hamiltonian and thus 



tr(E 



Efe-L H k 



0. Projections onto and orthogonal to 



Hff are most easily carried out by defining projectors 
P|| and P±. The actions of these projectors on the basis 
matrices <7j are given by 



P||<T, 



P±v x = 



2f2fcj Hj^ 



(14a) 



(14b) 



and the actions of P± on the remaining basis matrices 
are obtained by cyclic permutations. 

Secondly, SEk is likewise divided into two terms: Sem > 
commuting with the spin-orbit Hamiltonian and SEk±, 



orthogonal to it. It is helpful to think of Sem as the 
distribution of conserved spins. Ssk± can be thought of 
as the distribution of precessing spins. Equation (fT2|) is 
divided into separate equations for Sem and Seu±- 



dS 



Ek\\ 



dt 



P\\Jq (SEk) = ^Ek\\ 



J^" + T [Hk, SEk±] + P±Jo (SEk) = Efifcj 



(15a) 
(15b) 



The absence of the commutator [Hk, Sem] = in Eq. 
(|15ap indicates the absence of spin precession, while the 
commutator [Hk, SEh±] in Eq. (|15bj) represents spin pre- 
cession. In order to solve Eqs. (|15a|) and (|15b[) for ar- 
bitrary scattering, it is necessary to expand Ssk\\ and 
SEki. in the transition rate Wkk' , as 



S 



Ek\\ 



s 



Ek± 



= s 
= s 



(-1) 

Ek\\ 



s 



(0) 
Ek± 



(0) 
Ek\\ 

'EkJ 



q(1) 
D Ek\\ 

0(W, 



kk> 



<D{Wl k ,), (16a) 
(16b) 



This expansion is indeed an expansion in the parameter 
h/(Qk t p ), a fact that can be most clearly seen by exam- 
ining Eq. (|16bp and noting that the calculation of each 
term involves integration over time, which brings in a 
factor of l/fife, and the action of Jq. This expansion is 
therefore most suited to systems in the weak scattering 
regime. The expansion of Sem starts at order — 1, a fact 
which can be understood by inspecting Eq. (|15a[) . In the 
steady state the time derivative drops out, and the oper- 
ator Jo is first order in Wkk' , while the right-hand side 
is independent of Wkk'- As a result, the expansion of 
the solution must start at order —1. We examine next 
Eq. (|15b|) for SEk±- Since Hk is independent of Wkk', 
and the right hand side is also independent of Wkk' , the 
expansion of SEk± must start at order zero. 

Having divided the nonequilibrium correction to 
the spin density matrix into a part due to con- 
served spin and one due to precessing spin, we 
wish to determine the contributions these parts 
make to spin densities&^&^&^ i 10 ' 1 ^ 12 and spin 

currcn t 3 13 - 23 - 24 - 26 - 27 ■ 28 - 29.30. 31. 32. 33. 34. 35. 36. 37.38. 39. 40.41. 42.43. 44.45. 46.47 .4; 

in the steady state. The steady state density of spin 
component a is found by taking the trace tr(s CT SEk), 
where s a = (H/2)a CT . The steady-state spin cur- 
rent is found by taking the trace tr(j^f Seh), where 
Jf = hkjs^/m* + (l/4h)dn r7 /dk :j t. (The scalar term 
has zero expectation value.) Henceforth, for clarity 
and definiteness, integrals over wave vectors will be 
represented as two-dimensional. In the integrals below 
9' refers to the polar angle of k' . The extension to three 
dimensions is straightforward. 



B. Steady state for conserved spins 

We have shown that the expansion of S^kW begins at 
order —1, and we wish to find the first term in this ex- 
pansion. In the steady state, Eq. (|15a|) for the first term 
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in this expansion can be written as 

P \\Jo (S Ek ^) = %Ek\\- 
This equation can be recast as 



S 



(-1) 



where we have introduced 



J 0°Ek\\ ~ 



2-Kh 3 



2-kW> 



d6'W kk> 



de'w kk ,s { -v v 



(17) 



(18) 



(19a) 



To is the quantum lifetime of the charge carriers, i.e., the 
time between two consecutive scattering events. It differs 
from the momentum scattering time of Eqs. (fTTj) because, 
for nonisotropic scattering mechanisms, the information 
about the initial momentum is not lost after time to. 
Equation (TT5)) can be solved iteratively for any scattering 



(20) 



^Ek\\ ~ ^Bfc|| T + -f > ||^o(^'£;fc||) T 



P||J [P||J (S £fc ||)]r 3 



The equations for higher orders in W kk > are easily de- 
duced. However, the term of order —1 is by far the dom- 
inant one in the weak momentum scattering regime and 
is expected to be dominant over a wide range of strengths 
of the scattering potential. 

We examine more closely the nature of the steady state 
established for conserved spins. It is evident that this 
steady state involves no spin precession, and that the 
correction SE k \\ depends explicitly on the noncquilibrium 
shift in the Fermi surface and diverges in the ballistic 
regime, as To — > oo. In addition, it is important to 
note that scattering terms contain only the even func- 
tion W kk i . As a result, the correction S^kW does not give 
rise to a spin current. Inspection of Eq. (j2T)|) shows that 
integrals of the form 



d9 J° S Ek 



(21) 



contain an odd number of powers of k and are therefore 
zero. Consequently, in the absence of impurity spin-orbit 
interactions, the distribution of conserved spins can give 
no spin current. It can, however, give rise to a noncqui- 
librium spin density2*i^£&£&£ ' 1 1 1 1 1 1 2 , since integrals of 
the form 



dOS" S 



Ek\\ 



(22) 



C. Steady state for precessing spins 

The equations for the contributions to Sei^i. of orders 
zero and one in W kk i are 



dS 



(0) 

EkA 



dt 



X Y n hi a Ek^ 



} = ^Ek±+P±J (S E k\\), (23a) 



dS 



(i) 

Ek± 



dt 



[H k , S 



Ek±i 



(23b) 



To solve the equation for SeI± it is easiest to go into the 

interaction picture, obtain an expression for SeIj_, and 
then transform back to the Schrodinger picture. This 



(19b) procedure yields for S 



(o) 

Ek± 



s 



(o) _ 

Ek± ~ 



1 "fc x PW + p±J Oggfcji)] ■ g 

2 n k /h 



(24) 



where we have written Y<E k ±_ = (l/2)S_Efci_ ■ 17 an d 
SEk\\ = (1/2) Sem " °"- The result expressed by Eq. 
([2^)1 is valid for any elastic scattering. Since it does 
not depend explicitly on the form of the impurity po- 
tential, this term is usually regarded as intrinsic. Terms 
of higher order in W kk i are regarded as extrinsic because 
they depend explicitly on the form of the impurity po- 
tential. Nevertheless, it is evident from our work that, if 

+ P±Jq {Sem) vanishes, then vanishes and 

all the terms in Seic of higher order in W kk i also vanish. 

The steady state established for precessing spins is fi- 
nite in the clean limit. An argument similar to that given 
above for the distribution of conserved spin S^fcii shows 

that SeI± cannot lead to a nonequilibrium spin density 
(although, as will be shown below, higher-order terms in 
SEk± can contribute to the spin density). For, taking 
the expectation value of the spin operator, one arrives 
at integrals of the form J d9 s a Se\± , which involve odd 
numbers of powers of fc and are therefore zero. This term 
in the distribution of precessing spin does, however, give 
rise to nonzero spin currents, since integrals if the form 
J d6 J° SEk± contain an even numbers of powers of k 
and may be nonzero. Consequently, in the absence of 
spin-orbit coupling in the scattering potential, nonequi- 
librium spin currents arise from spin precession. 

In order to investigate terms of higher order in W kk i , it 
is easiest to examine a concrete case. Although the gen- 
eral conclusions apply to any elastic, spin-independent 
scattering, this analysis will be done in the next section 
in the context of short-range impurity scattering, where 
an exact solution is possible, which reveals an interesting 
physical picture. 



D. Short-range impurities 



contain an even number of powers of k and may be 
nonzero. 



An enlightening closed-form solution can be found for 
<Sgfc|| and SE k i_ for short-range impurities. In this case 



7 



Wkk' = W is a constant and J Ssfc = (S Ek ~ S E k)/r p , 
where the scattering time r p — H 3 /(mW) and the bar 
represents averaging over directions in k. Equation (|20p 
yields a closed- form solution for S Eh \\ 



S 



Ek\\ 



Pu (I - Pi 



(25) 



This solution enters Eq. (f24|) for 5 f ^ ± . The equations 
for the contributions to S E ki_ of higher orders in W can 
be easily determined 



'Ek± 



1 Hfc x |fE fc x [S Bfe± - J (SW||)]} 
2 



n 2 r 



(26a) 



o(2) _ 1 



i ri fc x {^ fc x {n k x [s Bfc± - j (s Ek \\)}}} ■ tr 



k 'p 



(26b) 



and so on. Since Y< E k± is orthogonal to the spin-orbit 
Hamiltonian HI, , we have ti(T, E k± Hf?) = Oj which tells 
us immediately that H E k±_ _L Jlfc. As a result, Jl^ x fi^ x 
^Ek± = —^Ek±, and 



1 EkA 



1 [SfifeJ 

2 



- Jo (5 



Efe|| 



5 



(2) 
£fc_L 



1 fife X [£_EfcJ_ — Jo (<S , £7fe|| )] • c 

2 fiH7^ ■ 



(27a) 



(27b) 



It is evident that the terms S^ff and give two 

separate geometric progressions. These progressions are 
easily summed to give for S E k± 



S 



EkJ_ — 



f2fc x (T, E kj_T p + P± S Ek \\) ■ cr r 2 



2fi(l 



(S_Efc_LT p + P± S Ek \\) 



1 



(28) 



Once again, if T. E k±_ t p + P± S E k\\ vanishes, then all the 
corrections to S E ki_ of order zero and higher also vanish. 

Let us analyze the two terms in S E ki_- By identifying 
terms in S E k± even and odd in fe, as was done above, 
it is evident that the first term on the RHS of Eq. I[28p 
leads to a spin current, but no spin density. The second 
term does not lead to a spin current, but it does produce 
a spin density. Nevertheless, this term tends to zero in 
the ballistic regime as well as in the strong momentum 
scattering regime, and it is not expected to be dominant. 

The closed-form solution found in this section shows 
that, in the absence of spin-orbit interactions in the im- 
purity potential, there is only one spin current, which in 
the weak momentum-scattering limit is independent of 
T p and in the strong momentum-scattering limit is oc t 2 . 

Bearing in mind that if Sglj_ vanishes all corrections of 
higher order also vanish, we conclude that, in the absence 
of spin-orbit interactions in the impurity potential, the 
distinction between intrinsic (disorder-independent) and 
extrinsic (disorder-dependent) spin currents is not useful. 



E. Steady state spin densities and currents 

Our work helps to understand the origins of nonequi- 
librium spin densities and spin currents in electric fields. 
The preceding sections illustrate the fact that nonequi- 
librium spin densities have two origins. The first, giving 
the dominant contribution, arises from the absence of 
precession. As charge carriers are accelerated, a fraction 
of their spin is conserved and produces a steady-state 
spin density^^£&£&£ ' 1 1 1 1 1 1 2 a process which has no 
analog in charge transport. Thus the dominant contri- 
bution to the noncquilibrium spin density in an electric 
field exists because in the course of spin precession a com- 
ponent of each individual spin is preserved. For an elec- 
tron with wave vector k, this spin component is parallel 
to f2fc. In equilibrium the average of these conserved 
components is zero. However, when an electric field is 
applied, the Fermi surface is shifted, and the average of 
the conserved spin components may be nonzero, as illus- 
trated in Fig. 1. This intuitive physical argument, to our 
knowledge absent to date from the literature, explains 
why the nonequilibrium spin density oc r^ 1 and requires 
scattering to balance the drift of the Fermi surface. It 
is interesting to note, also, that, although spin densities 
in electric fields require the presence of band structure 
spin-orbit interactions and therefore spin precession, the 
dominant contribution arises as a result of the absence 
of spin precession. 

An additional contribution, which vanishes in both the 
ballistic and the weak momentum scattering regimes, is 
associated with spin precession. This contribution arises 
from the term on the last line of Eq. (|2"8j) . The origin of 
this term can be understood by noting that, in the pres- 
ence of an electric field, the effective magnetic field about 
which a spin precesses changes slowly»22 (This is true in 
between scattering events.) This change is contained in 
the gauge- invariant crystal wave vector k = q — eEt/h. 
The fact that the effective magnetic field is changing 
slowly causes the spin to acquire a small component 
in the direction in which the effective magnetic field is 
changing. This component is proportional to the rate of 
change of the effective magnetic field and therefore, in our 
case, to E. It is associated with the 'flow' of ttk around 
the Fermi surface discussed by Shytov et alM This argu- 
ment explains why this term in the spin density vanishes 
in the clean limit as well as in the strong momentum 
scattering limit. For, in the absence of scattering, as k 
changes, the effective magnetic field will circle around the 
Fermi surface, and the component of the spin following it 
will average to zero. In the strong momentum scattering 
limit, on the other hand, the spin will not have time to 
acquire a component in the direction in which the effec- 
tive magnetic field is changing, due to the high frequency 
of scattering events. 

Furthermore, in the absence of spin-orbit interactions 
in the impurity potentials, spin currents are associated 
with displacement of spins. The relation between spin 
currents and spin precession was made explicit in the 
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ky 



E=0 




ky 



E>0 




k, 



FIG. 1: Effective field flk at the Fermi energy in the Rashba 
model " (a) without 63 (E = 0) and (b) with an external elec- 
tric field (E > 0). 



work of Sinova et al£2- 

We remark that both Sem and SEk± are invariant un- 
der time-reversal. As a result, the tensor characterizing 
the response of spin currents to electric fields is invariant 
under time reversal, whereas the tensor characterizing 
the response of spin densities to electric fields changes 
sign under time reversal, as expected in both casesi 24 i 64 



Interplay of conserved and precessing spin 
densities 



It is enlightening to compare the results obtained in 
the absence of scattering (i.e., the clean limit) with the 
results obtained when scattering is present. The aim is 
to obtain an understanding of the way scattering pro- 
cesses affect steady-state spin distributions in electric 
fields. This is done by comparing results obtained using 
the approach outlined in this paper with results obtained 
previously using Green's functions approaches, both for 
the case when scattering is not included and for the case 
in which scattering is taken into account. This process 
will aid us in identifying the information contained in 
vertex corrections to spin-related quantities in the frame- 
work of Green's functions-based theories of systems with 
spin-orbit interactions. The nature of this information is 
by no means obvious, and we will show that it has no 
analog in charge transport. 

In the absence of scattering, Eq. (|T5|) takes the form 



dS. 



Ek 



01 



[H k ,S, 



Ek\ 



eE dS Qk 
h ' dk ' 



(29) 



Comparison of Eqs. (fT2"|) and (f2"9"| shows that the driving 
term in the equation for the nonequilibrium spin distri- 
bution Ssk is renormalized by the term J s (n^fc), which 
accounts for spin-dependent scattering. 

In addition, Eq. (|23ap shows that scattering mixes the 
distributions of conserved and precessing spins. This is 
so because when one spin at wave vector k and precess- 
ing about Hfc is scattered to wave vector k' and pro- 
cesses about fifc/, its conserved component changes, a 
process which alters the distributions of conserved and 
precessing spin. Consequently, scattering processes in 



systems with spin-orbit interactions cause a renormaliza- 
tion of the driving term for the spin distribution, con- 
tained in Efifei, as well as scattering between the con- 
served and precessing spin distributions, described by 
P-i-Jb (^Ek\\ )■ Our analysis suggests that contributions 
due to these two processes are contained in vertex cor- 
rections to spin-dependent quantities found in Green's 
functions formalisms. 

In two dimensions, for Hamiltonians linear in wave vec- 
tor, we find that the renormalization term J s (riEk) does 
not contribute to the spin current for any elastic scat- 
tering. Therefore, the vertex correction to spin currents, 
found in other work^ i 34 ' 35 i 39 ' 40 i 46 represents only scatter- 
ing between conserved and precessing spin distributions. 
By noting that in Eq. (|28|) the correction to the source 
term in the equation for the precessing spin distribution 
has the form P± SEk\\ and thus depends on the steady 
state spin density, it becomes evident that the existence 
of a nonzero nonequilibrium spin density does affect the 
spin current. 

Furthermore, in three dimensions, for electrons in 
zincblende crystals, which are described by the fc 3 - 
Dresselhaus model, we find that, for short range impu- 
rities, the spin current obtained after inclusion of scat- 
tering is the same as when scattering is not included. It 
is known that the vertex correction to spin currents also 
vanishes 41 - in these systems. Noting that the steady-state 
spin density in zincblende crystals vanishes by symme- 
try, and more generally vanishes in any non-gyrotropic 
medium, 2 these observations reinforce our conclusion re- 
lating to the connection between vertex corrections to 
spin currents and the presence of a steady-state spin den- 
sity. We therefore expect vertex corrections to spin cur- 
rents to vanish in non-gyrotropic materials, in which no 
steady-state spin density is possible in an electric field. 2 



G. Comparison with previous work 

Our calculations for known cases give results in agree- 
ment with previous work. Firstly, our results agree with 
previous calculations of nonequilibrium spin densitiesi 4 ^ 
Furthermore, in two dimensions, for Hamiltonians lin- 
ear in wave vector, the spin current vanishes for short- 



range impurities, 



7.34.35.39.40.46 



as well as for small- angle 



scattering 42,^1 (We find that in fact it vanishes for any 
clastic scattering.—) For spin-orbit Hamiltonians charac- 
terized solely by one angular Fourier component N (Ref. 
l4lh the spin current oc N. 

Spin currents in systems in which band structure spin- 
orbit interactions are negligible were studied by En- 
gel et al.^ who demonstrated that spin currents in those 
circumstances are due to skew scattering. Skew scatter- 
ing appears as a term of third order in the scattering 
potential Ukk' (which must include spin-orbit coupling 
explicitly), whereas in this work we have restricted our 
discussion to terms of second order in Ukk', and we have 
not considered higher-order skew-scattering effects. 
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H. Observable effects 

Spin densities and spin currents excited by electric 
fields give rise to observable effects. In materials in which 
a nonequilibrium spin density is excited by an electric 
field, the presence of this spin density can be observed, 
for example, by means of magnetic circular dichroism 65 , 
which has long been used as a characterization tool for 
magnetic materials. Similarly, a spin current flowing 
transversely to the direction of the charge current (i.e., 
the spin-Hall effect) will give rise to a spin accumula- 
tion at the edge of the sample. A spin current flow- 
ing parallel to the direction of the charge current could 
be used as a means of spin injection from one semicon- 
ductor into another, as discussed in our recent worki^ 2 - 
Spin accumulation as a result of a spin-Hall current, or 
a spin density injected by means of a longitudinal spin 
current, can in turn be measured using magnetic circular 
dichroism^S, techniques developed recently. We note that 
alternative techniques can be used to observe nonequilib- 
rium spin densities^ ' 10 ! 11 ' 12 spin currents ; 19 ' 20 and edge 
spin accumulations . 16 ' 17 ! 18 

IV. SUMMARY 

We have demonstrated that, in the presence of band 
structure spin-orbit interactions, the steady state estab- 
lished for the carrier spin distribution contains two qual- 
itatively distinct contributions, corresponding to con- 



served and precessing spin. The distribution of conserved 
spin acquires a nonequilibrium correction that diverges 
in the ballistic regime. This correction is responsible for 
the establishment of the dominant nonequilibrium spin 
densities in electric fields. The distribution of precessing 
spin acquires a nonequilibrium correction that is finite in 
the ballistic regime. This correction is responsible for the 
establishment of nonequilibrium spin currents in electric 
fields and a small nonequilibrium spin polarization, which 
vanishes in the ballistic and strong momentum scattering 
regimes. We have demonstrated that, when spin-orbit in- 
teractions are present in the band structure and absent 
from the impurity potential, there is only one contribu- 
tion to the spin current, which appears independent of 
disorder in the ballistic regime but dependent on disorder 
in the strong momentum scattering regime. Moreover, 
we have also shown that scattering processes in systems 
with spin-orbit interactions give rise to a renormalization 
of the driving term in the equation for the spin distri- 
bution, as well as scattering between the conserved and 
precessing spin distributions, which sheds light on the 
nature of vertex corrections in these systems. 

The authors would like to acknowledge enlightening 
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ald, J. Sinova, D. L. Smith, J. Shi, H. A. Engel, and 
R. A Duine. The research at Argonne National Labo- 
ratory was supported by the US Department of Energy, 
Office of Science, Office of Basic Energy Sciences, under 
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APPENDIX A: MATRIX ELEMENTS OF A 
SCREENED COULOMB POTENTIAL 

In two dimensions, the matrix element Ukk' of a 
screened Coulomb potential between plane waves is given 
by 



Ukk' — — 



Ze 1 



1 



e V d=2 ^/\k-k'\ 2 + 1/L 2 .' 



(Al) 



where Z is the ionic charge, Vd=2 corresponds to the unit 
cell area, and L s is the screening length. The correspond- 
ing expression in three dimensions is 



Ukk' 



1 



Ze 2 



e V d=3 \k-k'\ 2 + l/L 2 ' 



(A2) 



where now Vd=3 is the unit cell volume. 



J 



APPENDIX B: INTEGRALS OVER TIME 

For a Hamiltonian given by Hk = ^ofel + (1/2) cr ■ fi*., the product of two time evolution operators e~ lHk,t e lHht 
can be written as 



e -iH k ,t e iH k t' _ e i{eo-e' )t' 



Clkt' ^'k^' r\ Qkt' ^'k^' ^ Qkt' ^fc^ 

cos cos — - i a ■ Ylk' cos sin — - Via- Ilk sin cos — — 

2 2 22 22 



(flk ■ fife' + ic ■ x ^fc) sin — — sin — ^ 



(Bl) 



with a similar expression holding for the Hermitian conjugate of this product. In the time integrals the trigonometric 
functions are expressed as complex exponentials, the integrals are evaluated, and the results are replaced by their 
principal parts, yielding 



dt > e i( S0 - s ' )s cos i^i cos -^ 



di J e i(e -e' )t' 



COS ■ 



■ sin - 



dt ' e ^o-e' )t' sin i^ cos "fe l 



dt > e i(zo-e' )t' sin ^i gfc^V 



\ Wk+ - e' k+ ) + 5(e k - - e' fe _) + 5(e k 



Ti Wk, 



6(e k - - e' k 



— [S(e k + - e' k+ ) + 5{e k+ - e' k _) - 5(e k - 



~k-) 



S (tk- - e'k 
<5( e fc+ - 4+) - S (tk- 



- 8(e k ~ 



,o 2 2 4 

Evaluation of all time integrals in this manner leads to Eq. ([7]) for the scattering term. 



k+) 



c k+) 



[S(fk+ - e'fc-) + #0fe- - e 'k+) - S ( e k- - e'k-) - S (ek+ - e k+ ) 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



